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ABSTRACT 



We report the discovery of HATS-2b, the second transiting extrasolar planet detected by the HATSouth survey. HATS- 
2b is moving on a circular orbit around a. V — 13.6 mag, K-type dwarf star (GSC 6665-00236), at a separation of 
0.0230 ± 0.0003 AU and with a period of 1.3541 days. The planetary parameters have been robustly determined using 
a simultaneous fit of the HATSouth, MPG/ESO 2.2m/GROND, Faulkes Telescope South/Spectral transit photometry 
and MPG/ESO 2.2m/FEROS, Euler 1.2m/CORALIE, AAT 3.9m/CYCLOPS radial-velocity measurements. HATS- 
2b has a mass of 1.37 ± 0.16 Mj, a radius of 1.14 ± 0.03 Rj and an equilibrium temperature of 1567 ± 30 K. The host 
star has a mass of 0.88 ± 0.04 Mq, radius of 0.89 ± 0.02 Rq and shows starspot activity. We characterized the stellar 
activity by analysing two photometric-follow-up transit light curves taken with the GROND instrument, both obtained 
simultaneously in four optical bands (covering the wavelength range of 3860 — 9520 A). The two light curves contain 
anomalies compatible with starspots on the photosphere of the parent star along the same transit chord. 



Key words, stars: planetary systems - stars: individual: (HATS-2, GSC 6665-00236) 
- techniques: spectroscopic, photometric 



■ stars: fundamental parameters 



1. Introduction 

^H 

;^ The first detection of a planet o rbiting a main-sequence star 
. ^ (51 Peg; Mayor & Queloz 1995 1 started a new era of astron- 
S^ omy and planetary sciences. In the years since, the focus on 
;_! exoplanetary discovery has steadily increased, resulting in 
Ci more than 850 planets being detected in 677 planetary sys- 
tem^ Statistical implications of the exoplanet discoveries, 
based on di fferent detection methods, have a.lso been pre- 
sented (e.g . [Mayor et al. 2011| Howard et al.||2012J Cassan 



et al. 112012 Fressin et al. 2013). Most of these planets have 



been detected by the transit and radial velocity (RV) tech- 
niques. The former detects the decrease in a parent-star's 
brightness due to the transit of a planet in front of it, while 
the latter measures the Doppler shift of parent-star light 



due to stellar motion around the star-planet barycenter. 
In the case of transiting extrasolar planets, the powerful 
combination of both methods permits a direct estimate of 
mass and radius of the planetary companion and therefore 
of the planetary average density and surface gravity. Such 
information is of fundamental importance in establishing 
the correct the oretical framework of planet formation a nd 
evolution (e.g. .Liu et al.||2011[ pordasini et~aLl|2012a|b ). 

Thanks to the effectiveness of ground- and s pace- 
based transit surveys lik e TrES ([Alonso ct al."200 4|), XO 
dMcCullough et a l.|2005t, HATNet (e. g. Bakos et al ,2012a 
Hartman et al. 201^ HATSouth jPenev et al|J2013p ' 



* Alfred P. Sloan Research Fellow 
*"* Packard Fellow 
^ exoplanet. eu, as at 2013, March 28 



W^ 



2012 



W ASP (e.g. [H cllicr e t al.||2012| |Smaliey et al-'PH^T'OES 
(Alsubai et al. 2011;|Bryan et al.||2012L KELT (ISiverd et 



aL||2012| ), COROT (e.g. [Rouan et al."2012'; •Patzold et aT 



) and Kepler ( Borucki et al. ]2011a_b , Batalha et a. 
), one third of the transiting exoplanets known today 
were detected in the past 2 years. In some cases, exten- 
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sive follow-up campaigns have been necessary to determine 
the correct physical properties of several planetary systems 



(e.g. [So uthworth et aL]|2011 Barros et al.|2011 Mancini et 
[al. 2013a|), or have been used to disco ver other planets b y 
measuring transit time variations (e.g. StefFen et al.||2013 |. 
With high-quality photometric observations it is also pos- 
sible to detect transit anomalies which a re connected with 
physical phenomena, such as star spots (|Pont et al. 2007 



Rabus et al |2009[ |pesert|201l][TVe gloan-R eed et al.|20T3| 



gravity dar kening (|Barnes||2009f Hzabo et al||201ip , stellar 
pulsations ( Collier Cameron et al. P!)10[ ), tidal distortion 
(Li et al. 2010 Leconte et al. 201 11), and the presence of 
additional bo dies (exomoons) (Kipping et al.|2009 Tusnski 
fc Valio|[20lT| ) 



In this paper we report the detection of HATS-2b, the 
second confirmed exoplanet found by the HATSouth tran- 
sit survey. HATSouth is the first global network of robotic 
wide-field telescopes, located at three sites in the Southern 
hemisphere: Las Campanas Observatory (Chile), Siding 
Spring Observatory (Australia) and H.E.S.S. (High Energy 
St ereoscopic Sy s tem) site (Namibia). We refer the reader 
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to 



Bakos et al. 



(|2013|, where the HATSouth instruments 
and operations are described in detail. HATS-2b is orbiting 
a K-type dwarf star and has characteristics similar to those 
of most hot-Jupiter detected so far. Photometric follow-up 
performed during two transits of this planet clearly reveals 
anomalies in the corresponding light curves, which are very 
likely related to the starspot activity of the parent star. 



2. Observations 

2.1. Photometry 

The star HATS-2 (GSC 6665-00236; V = 13.562 ± 0.016; 
J2000 a = 11'M6'"57".38, 5 = -22°21'46".77, proper mo- 
tion Ha = —45.8 i fc l.lmas/yr, us = ~ l-3 ± 1.5mas/yr; 
UCAC4 catalogue, Zacharias et al.|[2012 ) was identified as 
a potential exoplanet host based on photometry from all 
the instruments of the HATSouth facility (HS-1 to HS-6) 
between Jan 19 and Aug 10, 2010 (details are reported in 
Table fl]). The detection light curve is shown in Figure IT] 
This figure shows that the discovery data is of sufficient 
quality that it permits fitting a standard limb darkened 
planetary transit model. A detailed overview of the obser- 
vations, the data reduction and analysis is given in |Bakos| 
eFaLl(p013l). 
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Fig. 1. Top panel: Unbinned instrumental r-band light 
curve of HATS-2 folded with the period P = 1.354133 days 
resulting from the global fit described in Section [3] Lower 
panel: zoom in on the transit; the dark filled points show 
the light curve binned in phase using a bin size of 0.002. In 
both panels, the solid line shows the best-fit transit model 
(see Section 3.3). 



HATS-2 was afterwards photometrically foUowed-up 
three times by two different instruments at two differ- 
ent telescopes. On UT 2011, June 25, the mid-transit 
and the egress were observed with the "Spectral" imag- 
ing camera, mounted at the 2.0 m Faulkes Telescope South 
(FTS), at the Siding Spring Observatory (SSO). This par- 
tial transit permitted the refinement of the transit depth 
and ephemeris. The latter was particularly important for 
the subsequent follow-up observations performed with the 
MPGp]/ES0 2.2m telescope at the La Silla Observatory 
(LSOy Two full transits were covered on February 28 
and June 1, 2012, using GROND (Gamma-Ray Burst 
Optical/Near-Infrared Detector), which is an imaging 
camera capable of simultaneous photometric observations 
in four optical (ide ntical to Sloan g, r, i, z) passbands 
( Greiner et al.|[2008 ) . The main characteristics of the cam- 
eras and details of the data reduction are described in lPenevI 



et al. (2013). The GROND observations were performed 
with the telescope defocussed and using relatively long ex- 
posure times (80-90 s, 150-200 s cadence). This way min- 
imises noise sources (e.g. flat-fielding errors, atmospheric 
variation or scintillation, variation in seeing, bad tracking 
and Poisson noise) a nd delivers high-precisio n photome- 
try of transit events (Southworth et al. 2009). The light 
curves and their best-htting models are shown in Fig. 2j 
Distortions in the GROND light curves are clearly visible, 
which we ascribe to stellar activity. These patterns are anal- 
ysed in detail in Section 4. Table [T] gives an overview of all 
the photometric observations for HATS-2. 

2.2. Spectroscopy 

HATS-2 was spectroscopically followed-up between May 
2011 and April 2012 by five different instruments at five 
individual telescopes. The follow-up observations started in 
May 2011 with high signal to noise (S/N) medium resolu- 
tion (A/AA = 7000) reconnaissance observations performed 
at the ANU 2.3 m telescope located at SSO, with the im- 



agc slicing integral field spectrograph WiFeS (Dopita et al. 
,2007) . The results showed no RV variation with amplitude 
greater than 2kms~^; this excludes most false-positive sce- 
narios involving eclipsing binaries. Furthermore, an initial 
determination of the stellar atmospheric parameters was 
possible (Te//> = 4800 ± 300 if, logg^ = 4.4 ± 0.3), in- 
dicating that HATS-2 is a dwarf star. Within the same 
month, high precision RV follow-up observations started 
with the fibr e-fed echelle spectrograph CORALIE fQueloz| 
et al. 2000b I at the Swiss Leonard Euler 1.2 m telescope 
at LSO, followed by further high precision RV measure- 
ments obtaine d with the fibre-fed optic al echelle spectro- 
graph FEROS ( |Kaufer fc Pasquini|1998[ ) at the MPG/ESO 
2.2 m telescope at LSO. Using the s pectral synthesis code 



Max Planck Gesellschaft 



SME ('Spectroscopy Made Easy', Valenti & Piskunov 
1996) on the FEROS spectra, it was possible to deter- 



mine more accurate values for the stellar parameters (see 
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Table 1. Summary of photometric observations of HATS-2 



Facility 



Date(s) 



# of images Cadence (s) Filter 



Discovery 


















HS-1 (Chile) 






2010, 


Jan 24 - Aug 


09 


5913 


280 


Sloan r 


HS-2 (Chile) 






2010, 


Feb 11 - Aug 


10 


10195 


280 


Sloan r 


HS-3 (Namibia) 






2010, 


Feb 12 - Aug 


10 


1159 


280 


Sloan r 


HS-4 (Namibia) 






2010, 


Jan 26 - Aug 


10 


8405 


280 


Sloan r 


HS-5 (Austraha) 






2010, 


Jan 19 - Aug 


08 


640 


280 


Sloan r 


HS-6 (Austraha) 






2010, 


Aug 06 




8 


280 


Sloan r 


Follow-up 


















FTS/Spectral 






2011, 


June 25 




158 


60 


Sloan i 


MPG/ESO 2.2 m 


/ 


GROND 


2012, 


February 28 




69 


80 


Sloan g 


MPG/ESO 2.2 m 


/ 


GROND 


2012, 


February 28 




70 


80 


Sloan r 


MPG/ESO 2.2 m 


/ 


GROND 


2012, 


February 28 




69 


80 


Sloan i 


MPG/ESO 2.2 m 


/ 


GROND 


2012, 


February 28 




71 


80 


Sloan z 


MPG/ESO 2.2 m 


/ 


GROND 


2012, 


June 1 




99 


80 


Sloan g 


MPG/ESO 2.2 m 


/ 


GROND 


2012, 


June 1 




99 


80 


Sloan r 


MPG/ESO 2.2 m 


/ 


GROND 


2012, 


June 1 




99 


80 


Sloan i 


MPG/ESO 2.2 m 


/ 


GROND 


2012, 


June 1 




99 


80 


Sloan z 



i-band 



0.4 - 



-0.15 



2011 June25(FTS) 




2012 February 28 (GROND) 
f^-"--— 

2012 February 28 (GROND)- 



-0.1 



-0.05 0.05 

Time from transit center (days) 



0.1 



0.15 



Fig. 2. Phased transit hght curves of HATS-2. The dates 
and instruments used for each event are indicated. The hght 
curves are ordered according to the date and to the central 
wavelength of the filter used (Sloan g, r, i a nd z ). Our best 
fit is shown by the solid lines (see Section 3.3 1. Residuals 
from the fits are displayed at the bottom, in the same or- 
der as the top curves. The GROND datapoints affected by 
anomalies are marked with red empty squares and were not 
considered in estimating the final values of the planetary- 
system parameters (see Sect.[3|. 



Sect. 3.1 1. Further RV measurements were obtained with 
the CYCLOPS fibre-based integral field unit, feeding the 
cross-dispersed echelle spectrograph UCLES, mounted at 
the 3.9 m Anglo- Austrahan Telescope (AAT) at SSO, and 
with the fibre-fed echelle spectrograph FIES at the 2.5 m 
telesco pe at the Nordic Op tical Telescope in La Palma. We 
refer to Penev et al. (20131 for a more detailed description 
of the observations, the data reduction and the RV deter- 
mination methods for each individual instrument that we 
utilized. 

In total, 29 spectra were obtained, which are summa- 
rized in Table [2] Table [3] provides the high-precision RV 
and bisector span measurements. Figure [3] shows the com- 
bined high-precision RV measurements folded with the pe- 
riod of the transits. The error bars of the RV measurements 
include a component from astrophysical/instrumental jit- 
ter allowed to differ for the three instruments (Coralie: 
74.0 ms-\ FEROS: 44.0 ms^^, CYCLOPS: 193.0 ms"!). 

3. Analysis 

3.1. Stellar parameters 

As already mentioned in Sect. |2.2[ we estimated the stel- 
lar parameters, i.e. effective temperature Tf,ffi,, metallicity 
[Fe/H], surface gravity \ogg and projected rotational veloc- 
ity usini, by applying SME on the high-resolution FEROS 
spectra. SME determines stellar and atomic parameters by 
fitting spectra from model atmospheres to observed spec- 
tra and estimates the parameter errors using the quality of 
the fit, expressed by the reduced x^, as indicator. In case 
the S/N is not very high, or the spectrum is contaminated 
with telluric absorption features, cosmics or stellar emission 
lines, the reduced x^ does not always converge to unity, 
which leads to small errors for the stellar parameter values. 
To estimate of error bars, we used SME to determine the 
stellar parameters of each FEROS spectrum and calculated 
the weighted mean and corresponding scatter (weighted by 
the S/N of individual spectra). The results for the spectro- 
scopic stellar parameters of the SME analysis are listed in 
Table H 

By modeling the light curve alone it is possible to de- 
termine the stellar mean density, which is closely related 
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Table 2. Summary of spectroscopic observations of HATS-2. Observing mode: HPRV = high-precision RV measure- 
ments, RECON = reconnaissance observations. 



Tclcscopc/Instrunicnt 



ANU 2.3m/WiFcS 
Euler 1.2m/Coralie 
ESO/MPG 2.2m/FEROS 
ESO/MPG 2.2m/FEROS 
ESO/MPG 2.2m/FEROS 
ESO/MPG 2.2m/FEROS 
AAT 3.9m/CYCLOPS 
NOT 2.5m/FIES 



Date Range 



# of Observations Instrument resolution Observing mode 



2011, 


May 10-15 


2011, 


May 20-21 


2011, 


June 09-25 


2012, 


January 12 


2012, 


March 04-06 


2012, 


April 14-18 


2012, 


January 05-12 


2012, 


March 15 



7000 


RECON 


60000 


HPRV 


48000 


HPRV 


48000 


HPRV 


48000 


HPRV 


48000 


HPRV 


70000 


HPRV 


46000 


RECON 



Table 3. Relative RV and bisector span (BS) measurements of HATS-2 from various instruments used for high-precision 
RV measurements (c.f. Table [2]). 5 data points determined with FEROS are not listed here and were not used for further 
analysis due to high error bars caused by bad weather conditions. 



BJD 


Relative RV 


0"RV 


BS 


0"BS 


Phase 


Instrument 


(-2454000) 


(ms-i) 


(ms-1) 


(ms-i) 


(ms-1) 






1701.52346 


-90.37 


33.00 


-63.7 


79.5 


0.119 


Coralie 


1701.54622 


-238.37 


36.00 


117.7 


83.3 


0.135 


Coralie 


1702.52760 


168.63 


33.00 


31.3 


58.6 


0.860 


Coralie 


1702.55065 


169.63 


33.00 


-30.8 


74.2 


0.877 


Coralie 


1721.50300 


204.15 


71.76 


120.4 


93.8 


0.873 


FEROS 


1722.58300 


371.99 


100.77 


-146.9 


640.0 


0.671 


FEROS 


1723.44500 


-330.70 


73.66 


-211.0 


207.2 


0.307 


FEROS 


1736.46900 


164.63 


38.18 


149.3 


86.4 


0.925 


FEROS 


1737.53800 


261.69 


119.02 


-130.9 


287.5 


0.715 


FEROS 


1738.48600 


-47.00 


58.61 


-222.1 


146.9 


0.415 


FEROS 


1932.22448 


-33.63 


21.80 


1114.3 


43.2 


0.487 


CYCLOPS 


1933.21669 


-40.53 


51.20 


3464.1 


19.9 


0.219 


CYCLOPS 


1934.12774 


75.57 


65.58 


-4516.6 


168.1 


0.892 


CYCLOPS 


1938.81200 


-337.82 


139.13 


-436.4 


449.1 


0.351 


FEROS 


1939.16016 


47.07 


49.44 


9189.5 


1723.7 


0.608 


CYCLOPS 


1990.75600 


294.99 


103.71 


46.3 


201.6 


0.711 


FEROS 


1992.82100 


-335.90 


72.26 


-251.5 


208.5 


0.236 


FEROS 


2035.67400 


31.55 


77.56 


270.6 


223.1 


0.882 


FEROS 



Notes. The Coralie RV uncertainties listed here are known to be underestimated. Updated estimates are available, but we list 
here the values we used in the analysis. We note that in any case a jitter is included in the analysis to account for any additional 
scatter to that implied by the uncertainties (cf. Table [5|. 



to the normalized semimajor axis a/Ri, (Sect. 3.3) assum 



ing a circular orbit. Furthermore, adding RV measurements 
allows the determination of these parameters for elliptical 
orbits as well. 

To obtain the light curve model, quadratic limb- 
darke ning coefficien ts are needed, which were determined 
using Claret (20041 and the initially determined stellar 
spectroscopic para meters, ^ye use d the Yonsei-Yale stellar 
evolution models (Yi et al. 2001 hereafter YY) to deter- 



mine fundamental stellar parameters such as the mass, ra- 
dius, age and luminosity. The light curve based stellar mean 
density and spectroscopy based effective temperature, cou- 
pled with isochrone analysis, together permit a more accu- 
rate stellar surface gravity determination. The newly de- 
termined value for logg* = 4.50 ± 0.05 is consistent with 
the initial value of logg^ = 4.44 ± 0.12 thus we refrain 
from re-analyzing the spectra fixing the surface gravity to 
the revised value. The spectroscopic, photometric and de- 
rived stellar properties are listed in Table |4J whereas the 
adopted quadratic limb-darkening coefficients for the indi- 
vidual photometric filters are shown in Table [5] 

To illustrate the position of HATS-2 in the H-R dia- 
gram, we plotted the normalized semi-major axis a/ Ri, ver- 



sus effective temperature T^ffi,. Figure [4 shows the values 
for HATS-2 with their 1-a and 2-a confidence ellipsoids as 
well as YY-isochrones calculated for the determined metal- 
licity of [Fe/H]= 0.15 and interpolated to values between 1 
and 14 Gyr in 1 Gyr increments from our adopted model. 



3.2. Excluding blend scenarios 

To rule out the possibility that HATS-2 is actually a 
blended stellar system mimicking a transiting planet sys- 
tem we conduct a detailed modelin g of the light curves fol - 
lowing the procedure described in Hartman et al. (2011). 
Based on this analysis we can reject hierarchical triple star 
systems with greater than 4.5(7 confidence, and blends be- 
tween a foreground star and a background eclipsing binary 
with ^ 4cr confidence. Moreover, the only non-planetary 
blend scenarios which could plausibly fit the light curves 
(ones that cannot be rejected with greater than 5<7 confi- 
dence) are scenarios which would have easily been rejected 
by the spectroscopic observations (these would be obvi- 
ously double- lined systems, also yielding several kms~^ RV 
and/or BS variations) . We thus conclude that the observed 
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Fig. 3. Top panel: high-precision RV measurements for 
HATS-2 from CORALIE (dark fihed circles), FEROS (open 
triangles) and CYCLOPS (filled triangles) shown as a func- 
tion of orbital phase, together with our best-fit model. Zero 
phase corresponds to the time of mid-transit. The center-of- 
mass velocity has been subtracted. Second panel: velocity 
0-C residuals from the best fit. The error bars include a 
component from astrophysical/instrumental jitter allowed 
to differ for the three instruments. Third panel: bisector 
spans (BS), with the mean value subtracted. Note the dif- 
ferent vertical scales of the panels. 



oc 







HATS-2 


3.5 








;' / /'' /'' / / / 


4.0 




\ [ [ : ; : \ 


4.5 




\y\^>\\\SX\ 


5.0 




'!;'--y^>^SS<>:S$>:-~. 






--s^-^ ^^^-.'^~^r~~t'~-^"'-^>^v--^^l'^N>., 




--;;;;><:;-:-<2^^^^^^„ 


6.0 
6.5 




"~'-%;||i|||J|||5;,, 



5000 

Effective temperature [K] 



Yi et al. 



(2001) for the metal- 



Fig. 4. YY-isochrones from 
licity of [Fe/H]= 0.15. Isochrones are plotted for ages be- 
tween 1 and 14Gyr in steps of 1 Gyr (left to right). The 
ellipses mark the l-cr and 2-cr confidence ellipsoids for the 
determined values of a/Ri, and T^ffi,. 



Table 4. Stellar parameters for HATS-2 



Parameter 



Value Source 



Spectroscopic properties 

T.ff. (K) 

[Fe/H] ....^ 

vsini (kms ^) 

log g* (cgs) 

Photometric properties 

V (mag) 

B (mag) 

9 (mag) 

r (mag) 

i (mag) 

J (mag) 

H (mag) 

K (mag) 



5227 ±95 

0.15 ±0.05 

1.5 ±0.5 

4.44 ±0.12 



13.562 ± 
14.490 ± 
13.991 ± 
13.260 ± 
13.018 ± 
11.906 ± 
11.475 ± 
11.386 ± 



0.016 
0.031 
0.012 
0.020 
0.021 
0.024 
0.023 
0.023 



Derived properties 

M. (A/q) 

R* {Rg) 

log g* (cgs) 

L* (Lq) 

Mv (mag) 

Mk (mag) 

Age (Gyr) 

Distance (pc)"^ 



0.882 ±0.037 

0.898 ±0.019 

4.48 ±0.02 

0.54 ±0.06 

5.61 ±0.13 

3.64 ±0.07 

9.7 ±2.9 

360 ±11 



SME 
SME 
SME 
SME 



APASS^ 

APASS 

APASS 

APASS 

APASS 

2MASS2 

2MASS 

2MASS 



YY+a/i?.*+SME 
YY+a/iJ*+SME 
YY+a/iJ*+SME 
YY+a/R^+SME 
YY+a/fl*+SME 
YY+a/R^+SME 
YY+a/iJ*+SME 
YY+a/R^,+SME 



Notes. ^ The AAVSO (American Association of Variable Star 
Observers) Photometric All-Sky Survey 
^ Two Micron AU Sky Survey 
^ Av corrected 



transit is caused by a planetary companion orbiting HATS- 
2. 



3.3. Simultaneous analysis of photometry and radial velocity 



Following the procedure described in Bakos et al. (2010) 



the FTS and GROND photometric follow-up measurements 
(Tablell]) were simultaneously fitted with the high-precision 
RV measurements (Table l3|. The light curve parameters, 
RV parameters, and planetary parameters are listed in 
Tablell 

Table [5] also contains values for the radial velocity jitter 
for all three instruments used for high-precision RV mea- 
surements. They are added in quadrature to the RV results 
of the particular instrument. These values are determined 
such that x^ per degree of freedom equals unity for each in- 
strument when fitting a fiducial model. If x^ per degree of 
freedom is smaller than unity for that instrument, then no 
jitter would be added. The RV jitters are empirical numbers 
that are added to the measurements such that the actual 
scatter in the RV observations sets the posterior distribu- 
tions on parameters like the RV semi-amplitude. 

Allowing the orbital eccentricity to vary during the si- 
multaneous fit, we include the uncertainty for this value 
in the other physical parameters. We find that the observa- 
tions are consistent with a circular orbit (e = 0.071 ±0.049) 
and we therefore fix the eccentricity to zero for the rest 
of this analysis. Table [5] shows that the derived parame- 
ters obtained by including the distorted regions of the light 
curves are consistent with those derived with these regions 
excluded, indicating that the starspots themselves are not 
affecting the stellar or planet parameters in a significant 
way. We adopt the parameters obtained with the light curve 
distortions excluded in a fixed circular orbit. 
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Table 5. Orbital and planetary parameters for the HATS-2 system for different fitting scenarios: including the light 
curve distortions with free and fixed eccentricity e as well as excluding the light curve distortions with a fixed circular 
orbit. The last scenario was adopted for further analysis steps (parameters are highlighted in bold font). 



LC distortions included, 



LC distortions included, free e LC distortions excluded, e =: 



Parameter 

Light curve parameters 

P (days) 

Tc (BJD)" 

Ti4 (days)-^ 

Ti2 = T34 (days)-^ 

a/R^, 

i/R*" 

Rp/R* 

b = a cos i/R-j, 

i (deg) 

Limb-darkening coefficients'^ 

ag (linear term) 

bg (quadratic term) 

ar 

bj- 

li 

h 

az 

bz 

Radial velocity parameters 

K (ms-i) 

e cos uj 

esincj 

e 

^ (deg) 

RV jitter Coralie (ms~^) 

RV jitter FEROS (ms-i) 

RV jitter CYCLOPS (ms-^) .. 

Planetary parameters 

Mp (Mj) 

Rp (Rj) 

C{Mp,Rp)d' 

Pp (g/cm^) 

logflp (cgs) 

a (AU) 

Te, (K) 

0= 

(F) (lO^ergs-i cm-^)f 

Notes. 

" Tc'- Reference epoch of mid transit that minimizes the correlation with the orbital period. BJD is calculated from UTC. T14,: 
total transit duration, time between first to last contact; T12 — T34,: ingress/egress time, time between first and second, or third 
and fourth, contact. 

*" Reciprocal of the half duration of the transit used as a jump parameter in our MCMC analysis in place of a/ Rt,. It is related to 
a/Rt, by the expression C/^* = o,/R* ■ (27r(l + esincj))/(P\/T~ 

"^ Values for a quadratic law given separately for the Sloan g, r, 



1.354133 ±0.000001 

2455954.58576 ± 0.00008 

0.0859 ± 0.0004 

0.0104 ±0.0004 

t- C7+0.06 

26.48 ±0.06 
0.1317 ±0.0007 

(,214+0.061 

87.8 ±0.7 



0.7052 
0.1168 
0.4756 
0.2487 
0.3617 
0.2744 
0.2861 
0.2873 



272.2 ±30.5 
0.000 
0.000 
0.000 
0.000 
74.0 
44.0 
193.0 



1.369 ±0.158 

1.139 ±0.025 

0.11 

1.15 ±0.15 

3.42 ±0.05 

0.0230 ± 0.0003 

1567 ±30 

0.062 ±0.007 

1.36 ±0.11 



1.354133 ±0.000001 

2455951.87748 ±0.00009 

0.0859 ± 0.0004 

0.0107 ±0.0004 

5.65 ±0.32 

26.58 ±0.07 

0.1326 ±0.0008 

n 265+°''^3 

U.^DO_Q 075 

87.4 ±0.7 



0.7052 
0.1168 
0.4756 
0.2487 
0.3617 
0.2744 
0.2861 
0.2873 



278.7 ±33.0 

-0.033 ±0.052 

-0.023 ±0.060 

0.071 ± 0.049 

216 ± 77 

74.0 

44.0 

193.0 



1.397 ±0.171 

1.131 ±0.065 

-0.26 

1.20 ±0.28 

3.43 ±0.08 

0.0230 ± 0.0003 

1554 ±57 

0.064 ± 0.009 

1.32 ±0.20 



1.354133 ± 0.000001 

2455954.58576 ± 0.00009 

0.0862 ± 0.0004 

0.0109 ± 0.0005 

5.50 ± 0.09 

26.52 ± 0.07 

0.1335 ± 0.0010 



0.271 



-1-0.055 



-0.074 

87.2 ± 0.7 



0.7052 
0.1168 
0.4756 
0.2487 
0.3617 
0.2744 
0.2861 
0.2873 



268.9 ± 29.0 
0.000 
0.000 
0.000 
0.000 
74.0 
44.0 
193.0 



1.345 ±0.150 

1.168 ±0.030 

0.08 

1.05 ± 0.14 

3.39 ± 0.05 
0.0230 ± 0.0003 

1577 ± 31 
0.060 ± 0.007 

1.40 ± 0.11 



b^VT 



20101 



e^) (Bakos et al. 

^ _ _ _ „ , . * and z filters, i'hese values were adopted from the tabulations 

by [Claret (20041) according to the spectroscopic (SME) parameters listed in Table B] 

Correlation coefficient between the planetary mass Mp and radius Rp. ^_^ 

The Safronov number is given by 6 = (f/2)(Ksc/Vo,.6)^ = {a/Rp) ■ (Mp/M^) (see [Hansen fc Barman |2007 |. 
Incoming flux per unit surface area, averaged over the orbit. 



4. Starspot analysis 



4.1. Starspots and plages 



Fig.[5]shows the combined four-colour GROND light curves 
for the two HATS-2 transit events that were observed with 
this imaging instrument. The slight difference in the transit 
depth among the datasets is due to the different wavelength 
range covered by each filter. In particular, the g, r, i and z 
filters are sensitive to wavelength ranges of 3860 — 5340 A, 
5380 - 7060 A, 7160 - 8150 A, and 8260 - 9520 A, respec- 
tively. 



From an inspection of Fig. [51 it is easy to note several distor- 
tions in the light curves. Such anomalies cannot be removed 
by choosing different comparison stars for the differential 
photometry, and we interpret them as the consequence of 
the planet crossing irregularities on the stellar photosphere, 
i.e. starspots. It is well known that starspots are at a lower 
temperature than the rest of the photosphere. The fiux ra- 
tio should be therefore lower in the blue than the red. We 
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Fig. 5. Combined four-colour transit light curves of HATS-2 obtained with the GROND imaging system. Green dots 
are for the data taken in the g band, yellow upside down triangles for the r band, red squares for the i band, and purple 
triangles for the z band. To'p panel: transit observed on February 28, 2012. The bump observed just after the midtransit is 
interpreted as the covering of a "cold" starspot by the planet. Lower panel: transit observed on June 1, 2012. In addition 
to the bump occurred near the egress part of the light curve, a "hot" spot manifested in the g band, just before the 
starting of the covering of the starspot. 



thus expect to see stronger starspot features in the bluest 
bands. 

The data taken on February 28, 2012, are plotted in 
the top panel of Fig. [5] where the bump, which is clearly 
present just after midtransit in all four optical bands, is 
explained by a starspot covered by the planet. In particu- 
lar, considering the errorbars, the g, r, and i points in the 
starspot feature look as expected, whereas the feature in z 
is a bit peaked, especially the highest points at the peak of 
the bump at roughly BJD(TDB) 2455985.735. Before the 
starspot feature, it is also possible to note that the fluxes 
measured in the g and r bands are lower than those in 
the other two reddest bands, as if the planet were occult- 
ing a hotter zone of the stellar chromosphere. Actually, the 
most sensitive indicators of the chromospheric activity of a 
star in the visible spectrum are the emission lines of Call 
HA3968, KA3933, and Ha A6563, which in our case fall on 
the transmission wings of the g and r GROND passbands. 



Within the transit observed on June 1, 2012, whose data 
are plotted in the lower panel of Fig. [5j we detected another 
starspot, which occurred near the transit-egress zone of the 
light curve. Again, before the starspot feature, we note an- 
other "hotspot" in the g band, which has its peak at roughly 
BJD(TDB) 2456079.681. 

While these hotspot distortions could be caused by dif- 
ferential color extinction or other time-correlated errors 
(i.e. red noise) of atmospheric or astrophysical origin, there 
is also a physical interpretation that seems potentially plau- 
sible, i.e. that of a "plage". A plage is a chromospheric 
region typically located near active starspots, and usually 
forming before the starspots appear, and disappearing after 
the starspots vanish from a particular area. Accordingly, a 
plage occurs most often near a starspot region. As a matter 
of fact, in the GROND light curves, our plages are located 
just before each starspot. One can argue that the plage in 
the second transit is visible only in the g band, but this 
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can be explained by temperature fluctuations in the chro- 
mosphere, which causes a lack of ionized hydrogen, and by 
the fact that the Ca II lines are much stronger than the Ha 
line for a K-type star like IIATS-2. Another argument sup- 
porting this plage-starspot scenario is that, for these old 
stars, a solar-like relation between photospheric and chro- 
mospheric cycles is expected, the photospheric brightness 
varying in p hase with that of the chromosphere ( Lockwood 
eFaL][2007| . 



We note, however, that if these are plages they must 
be rather different from solar plages, which are essentially 
invisible in broad-band optical filters unless they are very 
close to the solar limb. Detecting a plage feature through 
a broad-band filter near the stellar center suggests a much 
larger temperature contrast and/or column density of chro- 
mospheric gas than in the solar case. 

4.2. Modelling transits and starspots 

We modelled the GROND transit light curves of HATS- 
2 with t he PRISIV^ and GEMC codes ( |Tregloan-Reed et 
ari|2013). The first code models a planetary transit over a 



spotted star, while the latter one is an optimisation algo- 
rithm for finding the global best fit and associated uncer- 
tainties. Using these codes, one can determine, besides the 
ratio of the radii Rp/ Ri,, the sum of the fractional radii, 
Tp + Vi, = {Rp + R*)/a, the limb darkening coefficients, the 
transit midpoint Tg, and the orbital inclination z, as well as 
the photometric parameters of the spots, i.e. the projected 
longitude and the latitude of their centres {6 and (j>), the 
spot size Tspot and the spot contrast Pspotj which is basi- 
cally the ratio of the surface brightness of the spot to that 
of the surrounding photosphere. Unfortunately, the current 
versions of PRISM and GEMC are set to fit only a single 
starspot (or hotspot), so we excluded the p-band dataset 
of the 2"*^ transit from the analysis, because it contains 
a hotspot with high contrast ratio between stellar photo- 
sphere and spot, which strongly interferes with the best- 
fitting model for the light curve. 

Given that the codes do not allow the datasets to be fit- 
ted simultaneously, we proceeded as follows. First, we mod- 
elled the seven datasets (1*'* transit: g, r, i, z; 2"*^ transit: r, 
i, z) of HATS-2 separately; this step allowed us to restrict 
the search space for each parameter. Then, we combined 
the four light curves of the first transit into a single dataset 
by taking the mean value at each point from the four bands 
at that point and we fitted the corresponding light curve; 
this second step was necessary to find a common value for 
To, i, 9 and (f). Finally, we fitted each light curve separately 
fixing the starspot position, the midtime of transit Tg and 
the system inclination to the values found in the previous 
combined fit. 

The light curves and their best-fitting models are shown 
in Fig. [61 while the derived photometric parameters for each 
light curve are reported in Table[6] together with the results 
of the MCMC error analysis for each solution. 

Comparing Table [5] with Table |6] we find that the fit- 
ted light curve parameters from the analysis described in 
Section |3.3| are consistent with the parameters that result 
from the GEMC+PRISM model, except for the inchna- 
tion which differs by more than 2a. As already discussed in 



Section [3. 3[ the joint-fit analysis was performed both con- 
sidering and without considering the points contaminated 
by the starspots, and the results are consistent with each 
other. So, our conclusion is that the spots themselves are 
not systematically affecting the stellar or planet parame- 
ters in a significant way; the differences in the inclination 
between GEMC and our joint fit are most likely due to 
differences in the modeling. 



4.3. Starspots discussion 

The final value for the starspots angular radii comes from 
the weighted mean of the results in each band and is 
''spot = 9.02° ± 0.30° for the starspot in the l*"* transit 
(spot #1) and rgpot — 19.16° ± 1.08° for the starspot in 
the 2"^^ transit (spot #2), with a reduced xl of 0.78 and 
0.49 respectively, indicating a good agreement between the 
various light curves in each of the two transits. We note 
that the error of the angular size of the spot #2 is greater 
than that of the spot #1. While it may be that spot 4f2 
is larger, we caution that its position near the limb of the 
star makes it size poorly constrained. 

The above numbers translate to radii of 13 145 ±552 km 
and 27923± 1729km, which are equivalent 0.05% and 0.2% 
of the stellar disk, respectively. Starspot sizes are in general 
estimated by doppler-ima ging reconstructions (i.e. Collier 



I Cameron 



1992 



to 22% of a stel 



Vogt et al.„1999J and their range is 0.1%o 
ar hemisphere, the inferior value bein g the 
detection limit of this technique (Strassmeier 2009). Our 



measurements are thus perfectly reasonable for a common 
starspot or for a starspot assembly, and in agreement with 
what has been found in other K-type stars. 

Starspots are also interesting in terms of how the con- 
trast changes with passband. In particular, we expect that 
moving from ultraviolet (UV) to infrared (IR) wavelengths 
the spot becomes brighter relative to the photosphere. 
Considering the starspot #2, its contrast decreases from 
r to z, even though this variation is inside the la er- 
ror (see Table [6]). Considering that IIATS-2 has an ef- 
fective temperature Teffi, = 5227 ± 95 K and model- 
ing both the photosphere and the starspot as black bod- 
ies jRabus et a;r| '2009' 'S anchis-Ojeda fc Winn]|2011| ), we 
used Eq. (1) of[Silva (200^ to estimate the temperature 
of the starspot ^2 in each band, obtaining the following 



values: T„ 



spot#2,r 



= 4916 ± 105 K, T 



spot#2,: 



4895 ± 121 



^ PRISM (Planetary Retrospective Integrated Star-spot 
Model). 



K and Tspot#2 z — 4856 ± 120 K. The weighted mean is 
T;pot#2 = 4891.5 ± 66.2 K. 

Unlike starspot #2, the spot contrasts for starspot #1 
are inconsistent with expectations. The spot is too bright 
in r relative to the other bandpasses, and too faint in z. 
If we estimate the starspot temperature in each band, we 
find rspot#i.g = 4345 ± 97 K, T3pot#i,r = 4604 ± 109 K, 
T;pot#i.. = 4318 ± 128 K and T,pot#i,. = 3595 ± 180 K. 
While the temperature in r is in agreement with those 
of g and i a,t 1 — 2a level, and slight differences could 
be explained by chromospheric contamination (filaments, 
spicules, etc.), the temperature in z seems physically inex- 
plicable. This effect is essentially caused by the z points 
at the peak of the starspot, at phase ~ 0.004 (see Fig. 6|, 
which are higher than the other points. However, one has 
also to consider that errorbars in this band are larger than 
those found in the other bands. This is due to the fact 
that, since the GROND system design does not permit 
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Table 6. Photometric parameters derived from the GEMC fitting of the GROND transit light curves. 



1^* transit 




Parameter Symbol 


g r I z 


Radius ratio Rp/R* 
Sum of fractional radii rp + r* 
Linear LD coefficient iti 
Quadratic LD coefficient U2 
Inclination (degrees) ° i 
Longitude of spot (degrees) '"''' 6 
Latitude of Spot (degrees) "■'' 
Spot angular radius (degrees) '^ rspot 
Spot contrast '^ Pspot 


0.1348 ±0.0011 0.1324 ±0.0011 0.13145 ± 0.00096 0.1352 ± 0.0010 

0.2204 ±0.0018 0.2232 ±0.0018 0.2149 ± 0.0016 0.2161 ± 0.0015 

0.749 ±0.060 0.593 ±0.051 0.352 ±0.057 0.298 ± 0.039 

0.171 ±0.018 0.296 ±0.043 0.218 ±0.028 0.144 ±0.024 

85.26 ±0.40 85.26 ±0.40 85.26 ±0.40 85.26 ±0.40 

5.78 ±0.45 5.78 ±0.45 5.78 ± 0.45 5.78 ±0.45 

76.52 ±1.94 76.52 ±1.94 76.52 ±1.94 76.52 ±1.94 

8.85 ±0.61 10.01 ±0.72 8.93 ± 0.65 8.72 ±0.50 

0.304 ±0.035 0.546 ±0.048 0.464 ± 0.052 0.251 ±0.52 


2'"'' transit 




Radius ratio Rp/R* 
Sum of fractional radii rp + r* 
Linear LD coefficient iti 
Quadratic LD coefficient M2 
Inclination (degrees) " i 
Longitude of spot (degrees) "'*' 6 
Latitude of Spot (degrees) "''' (f) 
Spot angular radius (degrees) '^ rspot 
Spot contrast '^ Pspot 


0.1356 ±0.0012 0.13411 ±0.00093 0.1307 ± 0.0011 
0.2108 ±0.0019 0.2022 ±0.0012 0.2039 ± 0.0017 
0.473 ±0.057 0.399 ±0.049 0.252 ± 0.044 
0.250 ±0.038 0.230 ±0.025 0.316 ± 0.049 
85.89 ±0.40 85.89 ±0.40 85.89 ±0.40 
35.26 ±1.20 35.26 ±1.20 35.26 ±1.20 
80.60 ±2.10 80.60 ±2.10 80.60 ±2.10 
20.14 ±1.49 17.79 ±2.17 18.28 ± 2.31 
0.753 ±0.046 0.780 ±0.054 0.789 ± 0.047 



Notes. 

° This is a common value and was found from the preceding fit of the combined data (see text). 

^ The longitude of the centre of the spot is defined to be 0° at the centre of the stellar disc and can vary from —90° to 90°. 

'^ The latitude of the centre of the spot is defined to be 0° at the north pole and 180° at the south pole. 

'^ Note that 90° degrees covers half of stellar surface. 

" Note that 1.0 equals to the surrounding photosphere. 



to chose different exposure times for each band, we are 
forced to optimize the observations for the r and i bands. 
Consequently, considering both the filter-transmission ef- 
ficiency and the color and the magnitude of HATS-2, the 
SNR in these two bands is better than that in z, for which 
we have larger uncertainty in the photometry. Taking these 
considerations into account, we estimated the final tem- 
perature of starspot #1 neglecting the z-band value, and 
obtaining rspot#i = 4425 ± 63 K. In Fig. [Tithe final values 
of the temperature contrast of the two starspots are com- 
pared wit h those of a sample of dwarf stars, which was re- 
ported by Berdyugina (2005). The derived contrast for the 
HATS-2 starspots is consistent wi th what is observe d for 
other stars. As already observed by Strassmeier ( 2009 ) , the 
temperature difference between photosphere and starspots 
can be not so different for stars of different spectral types. 
Moreover, in the case of long lifetime, the same starspot 
could bee n seen at quite different temperature (Kang & 
Wilson|1989j . It is then very difficult to find any clear cor- 
relation between starspot temperatures and spectral classes 
of stars. 

The achieved longitudes of the starspots are in agree- 
ment with a visual inspection of the light curves. The lati- 
tude of starspot #1, 76.52° ± 1.94°, matches well with that 
of starspot #2, 80.6° ± 2.1°, the difference being within la. 

Multiple planetary transits across the same spot com- 
plex can be used to constrain the alignment between the 
orbital axis of the planet a nd the spin axis of the star (e.g. 
Sanchis-Ojeda et al. (2011)). Unfortunately, from only two 
transits separated by 94 days we cannot tell whether or not 
the observed anomalies are due to the same comp lex. It is 
possible that they are. Following Solanki (2003) we esti- 
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Fig. 7. Spot temperature contrast with respect to the pho- 
tospheric temperature in several dwarf stars. Gray square is 
from spot ^ 2, black circle from spot #1, open circles from 
Berdyugina| (12005) . The name of the star and its spectral 
type is also reported for most of them. Note that some stars 
appear two times. 



mate a typical lifetime of ~ 130 days for spots of the size 
seen here. Moreover the rotation period of Prot = 31 ± 1 d 
inferred assuming they are the same spot is consistent with 
the value of Pi-ot = (30 ± 10 d) sin «^ estimated from the 
spectroscopically-determined sky-projected equatorial ro- 
tation speed. If they are the same spot complex, then the 
sky-projected spin-orbit alignment is A = 8° ± 8°, which is 
consistent with zero. We caution, however, that this value 
depends entirely on this assumption which could easily be 
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Fig. 6. Phased GROND light curves of HATS-2b compared to the best GEMC fits. The hght curves and the residuals 
are ordered according to the central wavelength of the filter used. Left panel: transit observed on February 28, 2012. Right 
panel: transit observed on June 1, 2012; due to to presence of the hotspot, the g band was not analysed with GEMC (see 
text). 



wrong. Continued photometric monitoring of HATS-2, or 
spectroscopic observations of the Rossiter-McLaughlin ef- 
fect, are necessary to measure the spin-orbit alignment of 
this system. 

5. Conclusions 

In this paper we have presented HATS-2b, the second 
planet discovered by the HATSouth survey. This survey 
is a global network of six identical telescopes l ocated at 
three different sites in the Southern hemisphere (Bakos et 



al. 2013) 



The parameters of the planetary system were 
estimated by an accurate joint fit of follow-up RV and 
photometric measurements. In particular, we found that 
HATS-2b has a mass of 1.369 ± 0.158^/,/, and a radius of 
1.139 ± 0.025i?j. To set this target in the context of other 
transit planet detections, we plotted 4 different types of cor- 
relation diagrams for the population of transiting planets 
(Fig. Is]) . We analysed the location of determined parame- 
ters for HATS-2b and its host star IIATS-2 in the following 
parameter spaces: planetary radius Rp vs. stellar radius 
i?*, planetary mass Mp vs. planetary equilibrium temper- 



ature Teq^p, planetary radius Rp vs. planetary equilibrium 
temperature T^q^p, planetary radius Rp vs. stellar effective 
temperature T^ffi,, and planetary radius Rp vs. planetary 
mass Mp . As illustrated in Fig. [8j the analysed parameter 
relations lie well within the global distribution of known 
exoplanets. 

Within each correlation diagram, at least one well char- 
acterized exoplanet can be found, whose parameters are 
consistent with those of IIATS-2(b) within the error bars. 
Looking on the correlation between planetary and stellar 
radius, HATS-2(b) is almost like HAT-P-37(b). Comparing 
the planetary equlibrium temperature and planetary mass, 
HATS-2b is similar to TrES-2b. The relation between plan- 
etary equilibrium temperature and planetary mass shows 
an agreement with WASP-32b, while the relation between 
stellar effective temperature and planetary radius points 
out that IIATS-2b agrees well with WASP-45b within the 
error bars. The focus on the planetary parameters radius 
and mass reveals a similarity to the transiting planet TrES- 
2. Comparing the atmospheres of exoplanets with similar 
physical parameters will be especially important to pursue 
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Fig. 8. Correlation diagrams for confirmed transit planets (exoplanet.eu, last updated Jan 10 2013). From top to bottom: 
planetary radius Rp vs. stellar radius Ri,, planetary mass Mp vs. planetary equilibrium temperature Teq.p, planetary 
radius Rp vs. planetary equilibrium temperature T^q^p, planetary radius Rp vs. stellar effective temperature T^ff^-,, and 
planetary radius Rp vs. planetary mass Mp. The position of the HATS-2 and HATS-2b parameters, respectively, are 
marked in red squares. Left panels give a global overview of the position of HATS-2 and HATS-2b in the population of 
known transit planets (black filled circles). Right panels illustrate a zoom in; green diamonds represent known exoplanets 
which fall within the errorbars of HATS-2(b) parameters. Short cuts: HP ^ HAT-P, W = WASP, T = TrES. 
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with e.g. the future ECHO space mission (Tinetti et al. 



2012 1. 

Very interesting is the detection of anomahes in the 
two muhi-band photometric-foUow-up hght curves obtained 
with the GROND imaging instrument. We recognize the 
anomahes as starspots covered by HATS-2b during the 
two transit events, a nd used PRISM and GEMC codes 
( Tregloan-Reed et al.j2013J to re-fit the transit light curves, 
measuring the parameters of the spots. Both the starspots 
appear to have associated hot-spots, which appeared in the 
1^* transit in the g and r bands, and only in the g band 
in the 2"^^ transit. These hotspots could be physically in- 
terpreted as chromospheric active regions known as plages^ 
which can be seen only in the GROND's bluest bands. We 
estimated the size and the temperature of the two starspots, 
finding values which are in agreement with those found in 
other G-K dwarf stars. 
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